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A proposed boundary condition accounting for shear layer effects within the Ffowcs 
Williams-Hawkings radiation module of the CDUCT-LaRC code is investigated. The 
development and numerical justification of the boundary condition formulation are re- 
viewed. An initial assessment of the effectiveness of the shear layer correction is con- 
ducted through comparison with experimental data. Preliminary results indicate that 
the correction provides physically meaningful modifications of the baseline predicted di- 
rectivity patterns. Trends of peak directivity steepening and shifting that appeared in 
predicted patterns were found to follow similar structures in measured data, particularly 
at higher radiation angles. 


Nomenclature 
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k, k x 

Mi, M2 
M d 
m 
it 
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R s 

it 

t 

Zl, Z2 


sound speed, m/s 

effective wave speed, m/s 

free space and axial wave numbers, 1/m 

mean flow Mach number in regions 1 and 2 

Mach number in the duct 

circumferential mode number 

unit normal to the shear layer 

annular duct radius 0 ou tcr — r inner ), m 

local shear layer radius, m 

fluid velocity, m/s 

modal trace velocity on shear layer surface, m/s 

acoustic impedance on shear layer surface adjacent to regions 1 and 2 


Symbols: 

p ambient density, kg/m 2 

co circular frequency, rad/s 
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I. Introduction 


The capability to efficiently predict acoustic radiation from aircraft engines with acceptable 
accuracy is a valuable tool in the development and assessment of low-noise concepts. The fan 
noise propagation and radiation code CDUCT-LaRC 1 includes an acoustic radiation module that 
is based on the Ffowcs Williams- Hawkings (FW-H) equation with penetrable data surfaces. 2 In 
the current implementation, the source surface is taken to be the inlet or exhaust planes for inlet or 
aft-fan cases, respectively. In many situations, this provides excellent results. However, there are 
limitations inherent in these source surface selections. The first issue can arise in both inlet and 
aft-fan cases when large off-axis radiation angles are of interest. In these situations, the shielding 
and edge effects of the engine cowl are not properly taken into account. This results in predictions 
that are essentially symmetric about the source surface. The second issue arises in aft-radiation 
cases and is associated with the refraction and possible amplification of sound waves as they prop- 
agate through the shear layer. This work focuses on the investigation of a proposed approach to 
account for shear layer effects within the framework of the CDUCT-LaRC code. The following 
discussion is broken into two main components: description of the model for shear layer refraction 
and application toward verification of the model through comparison with experimental data. 


II. Shear Layer Model 


The shear layer between the bypass flow and external stream can refract the sound waves 
radiated to the farfield. Radiation results can be improved by including this effect, as well as 
the reflection of the sound in the bypass region from the solid surface external to the bypass duct 
surrounding the core flow (see figure [T]). One way to achieve this is to extend the CDUCT-LaRC 
bypass duct propagation calculation, which is based on the parabolic approximation, beyond the 
exhaust plane and to apply a boundary condition (BC) similar to an acoustic liner BC on the shear 
layer. Although presented previously,^ a brief review of this boundary condition may be useful. 
The development is based on the satisfaction of two conditions on the shear layer: 

1 . The particle displacement must be continuous across the shear layer. 


2. The acoustic pressure must be continuous across the shear layer. 


The first result in the derivation involves the acoustic velocity normal to the shear layer. This 
follows a derivation of the liner BC by Myers 4 as presented by Farassat and Dunn. 5 The subscripts 
0 and 1 are used for the background and acoustic quantities. The particle displacement amplitude 
normal to the shear layer is also denoted by q 3 and the time dependence is harmonic with the factor 
e~ m . Following the notation of Farassat and Dunn 5 [eq. (3) with e = 1], one may write 

^ =-/(0<73 = fcfio -(no- Vno) + u{ -no ] 1 

= [q 3 no-(no-'Vuo) + ui-no] ( 1 ) 


where the superscripts 1 and 2 designate regions 1 and 2 in figure [I] it is the fluid velocity, and 
it is the unit normal to the shear layer. 

From equation [TJ one finds 
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and 


,,2 _ to+[nO-(nO-V«C))] 2 _! ^ 

Wl M 1^1^. 

i® + [no • (no-Vwo )] 1 

This result shows that if u\ n is found on the shear layer from the parabolic approximation method 
in region 1 , then u \ n on the external side of the shear layer may be calculated. 

Next, the derivation of the impedance, zi, on the side of the shear layer adjacent to region 2 is 
presented beginning with the model and assumptions^] First, the shear layer is taken to be infinitely 
thin, an assumption previously utilized in equation [3] Second, the barber pole pattern of the modes 
present in the flow region is assumed to persist in the flow region inside the shear layer. This 
barber pole pattern has a trace velocity on the shear layer surface in the direction normal to itself. 
Figure [ 2 ] shows this trace velocity, V s , and the method by which it is calculated. This figure also 
shows a local tangent plane to the shear layer surface which is cylindrical. Assuming the acoustic 
pressure in a mode to be proportional to exp [— i (to t — m0 — k x x)] where m is the circumferential 
mode number, x is the axial distance, and k x is the axial wave number, it follows that 


6 = 



(4) 
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( 6 ) 


V s = 




(7) 


Here, R s is the local radius of the shear layer (assumed almost constant), k = to/ci, and ci is the 
local speed of sound in region 1 of the shear layer. 

From Morse and and Ingard 6 (equation 11.1.21), the impedance, z.2, is 


Z2 = 


P2C2 

(1 T M 2 cos (])) sin(]) 


(8) 


where M 2 = «2 / C 2 is the Mach number based on the background flow in region 2 and <|) is the 
preferred direction of propagation of plane waves in region 2. This angle is found purely by 
kinematic considerations. Figure [3 shows the geometry in the local plane containing and the 
local normal to the shear layer which is along axis 2. Let = (sin<|), cos(|)) be the direction normal 
to the wavefront. The effective wave speed in region 2, c^ , is 


~cl = C2~t + «2 COS\(r (9) 

where cos\|/ = x/V s is shown in figure |2] The projection of cj on the shear layer is equal to the 
trace speed V s . This gives 

V s = U 2 cost)/ + C 2 sin <}) 

a Note that until the end of this section, we work with the zeroth order quantities (background flow, time independent 
quantities). For simplicity of notation, we denote regions 1 and 2 by subscripts. 
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from which it is seen that 


. V s -u 2 cost)/ 

sintp = 

C 2 

Using this in equation [8] leads to 


P2C2 (Vs~U 2 ) 

V s \Z(V s -u 2 ) 2 -4 


The final step in the formulation is that 
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( 10 ) 


( 11 ) 


( 12 ) 


and equation |3] provides an expression for u\ n /u\ n . 

III. Implementation 

Implementation of the impedance boundary condition rests on the calculation of two main 
quantities. First, the ratio of particle velocity external to the shear layer (region 2) to the particle 
velocity internal to the shear layer (region 1), u\ n /u\ n must be known (see eq. [3]). Determina- 
tion of this ratio requires knowledge of the background flow velocity vectors and the unit normal 
to the shear layer surface. The background flow velocity vectors and location of the shear layer 
are available from the mean flow solution supplied for the propagation calculations. The Non- 
Uniform Rational B-Splines (NURBS) techniques utilized within the grid generation and redistri- 
bution functionality of CDUCT-LaRC provide the ability to define a general shear layer surface 
and therefore, the unit normal to the shear layer. 

Secondly, the impedance on the external side of the shear layer, zi, must also be specified. As 
seen in equation [TT] calculation of zn requires knowledge of the local radius of the shear layer, R s , 
the circumferential mode number, m, and the axial wave number, k x . The possible issues involved 
in obtaining these quantities in general has been discussed previously. 3 For the cases considered 
in this work, the possible complications are avoided by using the value of 22 obtained for the lowest 
radial order of the highest order circumferential mode. 

With the required quantities in hand, the shear layer model may incorporated into the propaga- 
tion module of CDUCT-LaRC. The shear layer surface defines the outer surface of a ’pseudo-duct’ 
that extends from the exhaust plane of the bypass duct. An inner ’pseudo-duct’ surface must also 
be defined, which in a typical engine application is at least partially defined by the outer surface 
of the core cowl. The source plane of the ’pseudo-duct’ is defined by the exit surface of the final 
block used in the internal propagation calculation. The acoustic potential values on this exit plane 
are also used directly to specify the generic source and an additional propagation calculation is 
subsequently performed within the ’pseudo-duct’. The exposed surfaces of the ’pseudo-duct’ are 
then used as FW-H source surfaces in the radiation calculations. In order to illustrate the imple- 
mentation further, and provide some verification of the approach, some sample calculations are 
presented. 
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IV. Application 


The verification cases considered are the by-product of previous calculations performed in con- 
junction with colleagues at NASA Glenn Research Center. In this previous work, CDUCT-LaRC 
predictions were made for the NASA Glenn Active Noise Control Fan (ANCF) rig and compared 
with results obtained using the Eversman finite element code. 7 The ANCF rig is a low speed fan 
testbed designed for noise reduction concept testing and code verification/development. 8 In addi- 
tion to collection of farfield data, the ANCF allows for the insertion of a rotating microphone rake 9 
to measure the interaction modes of the ducted fan. Therefore, for selected test configurations, in- 
duct rotating rake measurements and farfield acoustics measurements may be obtained. A general 
schematic of the axi-symmetric rig is presented in figure |4| 

Although not a focus of this work, some points may be made concerning the preliminary com- 
parisons between CDUCT-LaRC and the Eversman finite element code. First, as a result of initial 
problem formulation, the capability to directly input rotating rake data was added to CDUCT- 
LaRC. Secondly, to illustrate the effects of including the proposed shear layer correction, compar- 
ison of the CDUCT-LaRC and Eversman inlet and aft predictions for the (2,0) mode at one BPF 
(533 Hz) are provided in figure [5} Aft baseline and shear layer corrected CDUCT-LaRC predic- 
tions are included. Note that for these aft results, inclusion of the shear layer correction increases 
the peak directivity level and shifts its location forward. Overall, the initial comparisons were en- 
couraging and provided motivation to pursue further investigation of this approach to accounting 
for shear layer effects. Having introduced the ANCF configuration and some of the initial effort, 
comparisons with measured data for aft radiation are discussed. 

A range of source frequencies (1,2, and 3 BPF) and flow speeds (1700 through 2000 RPM) 
were considered. The corresponding rotating rake data for each case was used to specify the modal 
sources. The expected Tyler-Sofrin modes (based on 16 blades and 14 stators), with amplitude and 
phase for 1, 2, and 3 BPF at each fan speed, are listed in tables [T] through [4] of the Appendix. 
Note that the source information corresponding to the various BPF values are separated within 
the tables and the duct Mach number, Mj, is supplied as well. The rotating rake data is supplied 
at an axial location of x = 61.6 in. The proximity to the exhaust plane (x = 64.6 in.) leads to 
the use of a short duct section for the internal propagation calculations. Subsequently, to utilize 
the shear layer boundary condition, an external ’pseudo-duct’ must be attached to the internal duct 
segment. Propagation calculations are then performed within the ’pseudo-duct’ and the solution on 
the exposed surfaces used as the source surface for the FW-H radiation calculations. The farfield 
directivity obtained in this manner may depend on the axial extent of the ’pseudo-duct’. Therefore, 
the calculations were made for several ’pseudo-duct’ lengths defined in terms of multiples of the 
annular duct radius, R, at the exhaust plane. Specifically, lengths equivalent to 1.0, 2.5, 5.0, and 
10.0 multiples of the duct radius were investigated. As an example of the computational meshes 
used, the internal and ’pseudo-duct’ meshes for the 2.5 R case are shown in figure[6]and represented 
with white and green wireframe, respectively. Note that, although the geometry is axi-symmetric, 
the calculations are performed over the fully three-dimensional geometry for verification purposes. 

In this figure, it can also be seen that the shear layer is assumed to form a straight cylindrical 
surface. As such, the Mach number in the duct (internal and pseudo) is taken to be uniform and the 
external Mach number is set to zero. The aforementioned flow and source specifications were used 
to perform internal and ’pseudo-duct’ propagation calculations. For the internal propagation calcu- 
lation, the grid size was 21X89X9 in the radial, circumferential, and axial directions, respectively. 
All ’pseudo-ducts’ used the same radial and circumferential dimension as the internal block. The 
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axial grid sizes and typical wall-clock times for complete propagation and radiation predictions for 
each of the cases are provided in table [5j In this table, the EP case identifier refers to the baseline 
radiation calculation using the exhaust plane of the bypass duct as the FW-H source surface. Sub- 
sequent cases involving the shear layer boundary condition are referred to as SL:XX where XX 
refers to the length of the ’pseudo-duct’ in multiples of the annular duct radius. Farfield data for the 
ANCF is available at 15 microphone locations evenly spaced along a 15 foot arc centered on the 
exhaust plane. Figure [7] shows the placement of the observer arc relative to the ANCF geometry. 

As mentioned previously, predictions were performed for the first three blade passage frequen- 
cies for corrected fan speeds of 1700 to 2000 RPM. The results for each case are presented in 
figures [8] through [19} Note that, for convenience, the remaining figures are placed at the end of the 
document. Looking first at the 1 BPF results, it can be seen (figures[8j[TTJ[T4} and[l7]) that the effect 
of the shear layer boundary condition is to increase and narrow the peak directivity, as well as shift 
it forward. This trend compares favorably with the measured data and the predictions match quite 
well at 2000 RPM. The more complex directivity patterns of the 2 BPF cases ((9} [12} [15j and |T8j) 
result in greater disparity between the predicted and measured farfield values. However, the shear 
layer correction tends to modify the fairly smooth baseline result and provide reasonable detection 
of peaks that begin to appear. In particular, the results at 1700 and 1800 RPM (figures |9j[T2|) show 
peaks in the modified directivity patterns that are similar to the measured data, particularly at the 
higher radiation angles. Finally, the comparisons between predicted and measured data at 3 BPF 
(figures [10} [T3j [16} and [T9]) are similar to those for the 2 BPF cases. The shear layer correction 
produces directivity patterns that are more directional relative to the baseline prediction. Addi- 
tionally, the locations of the peaks appear to follow trends in the measured data, particularly at the 
higher radiation angles. Generally, the results are encouraging in that the shear layer correction 
produces modifications to the baseline predictions that more closely follow trends in the measured 
data. In several instances, the predicted directivity patterns show main peaks that are very much 
in line with the measured data in both angle and amplitude. Some of the discrepancies at forward 
radiation angles may be due to diffraction at the trailing edge of the nacelle. In the 2 BPF, 2000 
RPM case however, some of the difficulties may be partially related to the source description. As 
seen in table [4j the rotating rake data included information for the in = 4 an -8 circumferential 
modes. However, from a Tyler-Sofrin analysis, a circumferential mode order of -10 is possible 
due to interaction at 2 BPF. At 2000 RPM, the source frequency is such that the (-10,1) mode is 
just above cut-off. To gain some qualitative understanding of the effects of including this mode, a 
sample calculation was performed for the 2 BPF, 2000 RPM case in which the (-10,1) mode was 
added to the original source. The amplitude of the mode was set to 2.0 N/m 2 and the phase was set 
to zero. For illustrative purposes, the predicted pattern was offset by 6 dB and also plotted (SLA: 
10R). From comparison of figures [18} and |20} it can be seen that in this case that the inclusion of 
the (-10,1) mode greatly improves the comparison with measured data. For the lower fan speeds 
however, the (-10,1) is cut-off and the initial comment concerning trailing edge effects may be 
appropriate. Continued comparison with other codes, such as the Eversman finite element code, 
will provide further information on this issue. Additionally, further understanding may be gained 
by comparison with results presented by Mani 10 and MuntP^ 
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V. Concluding Remarks 


In this work, the motivation for, and description of, a proposed boundary condition for the FW- 
H acoustic radiation module of CDUCT-LaRC have been presented. Previous CDUCT-LaRC pre- 
dictions for the ANCF rig have compared favorably with results obtained by colleagues at NASA 
Glenn using the Eversman finite element code. This led to further study of the effectiveness of a 
shear layer correction developed for inclusion in the FW-H radiation module of CDUCT- LaRC. 
Preliminary hardwall ANCF predictions were generated, using rotating rake data as source input, 
and compared with measured farfield data. Results indicate that the correction provides physically 
meaningful modifications of the baseline predicted directivity patterns. Some discrepancies at for- 
ward radiation angles were improved by inclusion of additional source modes, providing further 
illustration of the known importance of source specification. Some of the other discrepancies at 
forward radiation angles may be the results of edge effects encountered at the trailing edge of 
the nacelle. Further comparison with results from the Eversman finite element code, as well as 
work by Mani and Munt, will provide further insight into this matter. Overall, predicted trends of 
peak directivity steepening and shifting were found to follow similar structures in measured data, 
particularly at higher radiation angles. 
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Figure 1. Boundary Regions for the Aft-Fan Geometry 
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Figure 3. Effective Wave Speed Description 
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Figure 4. Schematic of NASA Glenn ANCF rig 



Figure 5. Initial Comparison of CDUCT-LaRC and Eversman Predictions [BPF = 533 Hz, (2,0) mode] 
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Figure 6. Internal and ’Pseudo-Duct’ Computational Meshes (2.57?) 



Figure 7. Observer Arc Location Relative to the ANCF geometry 
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Figure 8. Aft ANCF prediction: 1700 RPM, 1BPF (453 Hz) 



Figure 9. Aft ANCF prediction: 1700 RPM, 2BPF (906 Hz) 
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Figure 10. Aft ANCF prediction: 1700 RPM, 3BPF (1360 Hz) 



Figure 11. Aft ANCF prediction: 1800 RPM, 1BPF (480 Hz) 
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Figure 12. Aft ANCF prediction: 1800 RPM, 2BPF (960 Hz) 



Figure 13. Aft ANCF prediction: 1800 RPM, 3BPF (1440 Hz) 
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Figure 14. Aft ANCF prediction: 1900 RPM, 1BPF (507Hz) 



Figure 15. Aft ANCF prediction: 1900 RPM, 2BPF (1013 Hz) 
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Figure 16. Aft ANCF prediction: 1900 RPM, 3BPF (1520 Hz) 



Figure 17. Aft ANCF prediction: 2000 RPM, 1BPF (533Hz) 
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Figure 18. Aft ANCF prediction: 2000 RPM, 2BPF (1066 Hz) 



Figure 19. Aft ANCF prediction: 2000 RPM, 3BPF (1600 Hz) 


16 of |T9] 


American Institute of Aeronautics and Astronautics 



Figure 20. Aft ANCF predictions: (-10, 1) mode addition [2BPF = 1066 Hz] 


Tables 


Mode 

Amplitude [N/m 2 ] 

Phase [rad] 

(2,1) 

5.12 

3.11 

(2,2) 

0.56 

1.74 

(4,1) 

4.53 

2.41 

(4,2) 

1.07 

-2.67 

(4,3) 

0.47 

-1.01 

(6,1) 

1.27 

1.68 

(6,2) 

0.70 

1.58 

(6,3) 

0.30 

1.46 

(6,4) 

0.28 

1.45 

(-8,1) 

1.68 

1.73 

(-8,2) 

0.83 

-0.75 

(-8,3) 

0.18 

-2.37 


Table 1. Source Information for 1700 RPM (Mj = 0.144) 
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Mode 

Amplitude [N/m 2 ] 

Phase [rad] 

(2,1) 

5.16 

2.08 

(2,2) 

0.91 

0.61 

(4,1) 

6.16 

-0.92 

(4,2) 

2.38 

0.15 

(4,3) 

0.46 

2.35 

(6,1) 

7.40 

0.86 

(6,2) 

0.47 

-0.27 

(6,3) 

1.51 

0.92 

(6,4) 

0.42 

-2.34 

(-8,1) 

3.13 

1.47 

(-8,2) 

0.62 

-3.04 

(-8,3) 

0.12 

-1.32 


Table 2. Source Information for 1800 RPM (Mj = 0.154) 


Mode 

Amplitude [N/m 2 ] 

Phase [deg] 

(2,1) 

7.82 

2.14 

(2,2) 

1.44 

1.70 

(4,1) 

2.59 

-0.62 

(4,2) 

1.27 

1.75 

(4,3) 

1.67 

2.45 

(6,1) 

7.70 

-0.64 

(6,2) 

0.93 

-2.81 

(6,3) 

0.96 

3.04 

(6,4) 

0.45 

1.80 

(-8,1) 

2.11 

0.44 

(-8,2) 

1.19 

-2.53 

(-8,3) 

1.33 

0.34 


Table 3. Source Information for 1900 RPM (Mj = 0.164) 
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Mode 

Amplitude [N/m 2 ] 

Phase [deg] 

(2,1) 

8.14 

0.04 

(2,2) 

2.01 

-0.80 

(4,1) 

1.91 

2.73 

(4,2) 

1.19 

2.07 

(4,3) 

0.28 

0.67 

(6,1) 

1.97 

-0.99 

(6,2) 

1.48 

2.10 

(6,3) 

0.49 

0.88 

(6,4) 

0.45 

1.81 

(-8,1) 

3.29 

-1.08 

(-8,2) 

2.13 

0.92 

(-8,3) 

2.27 

-0.89 

(-8,4) 

0.43 

- 0.72 


Table 4. Source Information for 2000 RPM (M ( j = 0.175) 


Case 

Pseudo-Duct Length 

Axial Grid Dimension 

Total Wall-Clock Time [s] 

EP 

0 

0 

3 

SL: 1.0R 

R 

33 

13 

SL: 2.5R 

2.5R 

81 

27 

SL: 5.0R 

5. OR 

161 

51 

SL: 10.0R 

10.0R 

321 

103 


Table 5. Grid Sizes and Complete Calculation Timings for Various ’Pseudo-Ducts’ 
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